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Abstract 

We present here fully optimized two-dimensional pupil apodizations for which 
no specific geometric constraints are put on the pupil plane apodization, apart from 
the shape of the aperture itself. Masks for circular and segmented apertures are 
displayed, with and without central obstruction and spiders. Examples of optimal 
masks are shown for Subaru, SPICA and JWST. Several high-contrast regions are 
considered with different sizes, positions, shapes and contrasts. It is interesting to 
note that all the masks that result from these optimizations tend to have a binary 
transmission profile. 

1 Introduction 

As the population of known exoplanets has grown over the last decade, interest has been 
steadily increasing in methods and missions to observe them directly. These include 
both ground and space based systems. In fact, the 2010 decadal survey for astronomy 
and astrophysics 1 1 1 placed a high priority on developing both the technology and mis- 
sion concepts for a high-contrast planet imager that would be ready for the next decade. 
The goal is to be able to image Earthlike planets in the habitable zones of their parent 
stars with a telescope from 4 to 8 m in diameter. Many approaches to achieving the 
high contrast needed have been proposed, including many types of coronagraphs, pupil 
interferometers, and external occulters. In this paper, we present new shaped-pupil de- 
signs that increase the system throughput and can be used with on-axis and segmented 
telescopes, introducing the potential for dramatically simplifying the engineering of a 
future ground or space telescope. 

Pupil apodization creates a high-contrast point-spread function by varying the am- 
plitude transmission in the entrance pupil of the telescope (see, e.g., |4|). In fact, 
Slepian [12] found the best possible apodization in the sense that it maximally concen- 
trates the light in the image plane toward the center leaving very high contrast outside 
the inner working angle. Unfortunately, manufacturing aperture masks with variable 
transmission over broad bands with sufficient accuracy has proved formidable. Instead, 
we have focused on binary masks that reshape the entrance pupil by either allowing or 
blocking transmission |7|. 

Our approach to designing shaped pupils has been to maximize the transmission 
of the pupil subject to constraints on the field distribution in the image plane. Be- 
cause a linear operator (the Fourier transform) is used to switch from one plane to the 



1 



other, these optimization problems are Hnear programs, and can thus be solved very 
efficiently. However, so far, full optimizations (i.e. without any constraints put on the 
pupil except for the initial transmission of the aperture) have only been done for one- 
dimensional geometries. The results of these optimizations include barcode masks, 
checkerboard masks and concentric ring masks ||9] [TOl \14\ . In the last two cases, a 
square geometry and a circular symmetric geometry, respectively, were assumed so 
that these ID problems resulted in 2D masks. 

An important observation is that each of these ID optimizations tended toward bi- 
nary distributions of amplitude in the pupil plane; their binary aspect is not a constraint 
but rather is a result of the optimization. The outer working angle (OWA) of these 
masks govern (indirectly) the number of the blocking and transmitting regions, and in 
the case of a large OWA the width of these regions becomes smaller than the width 
associated with one element of resolution. In this case the apodization profile is no 
longer binary but rather smooth. 

In this paper, we describe a new optimization approach that results directly in 2D 
shaped pupils with no simplifying geometric assumptions. The idea, first described in 
iTlSJ . exploits the main idea behind the fast Fourier transform (see also |17|) thereby 
allowing us to design masks for arbitrary pupil geometries, including those with central 
obstructions, spiders, or segmented telescopes. It also allows us to tailor the region of 
high contrast in the image plane to any desired shape or size. Thus, for instance, we can 
design for dark holes that match those created by the accompanying wavefront control 
system. This ability often results in systems with higher transmission, smaller inner 
working angles (IWA) or higher contrast. 

The paper is organized as follows. In Section|2j we describe the formalism for the 
optimizations. In Section[3] we present the result of the optimization for a circular un- 
obstructed aperture. We use this example to show how the sampling of the pupil plane 
affects the quality of the optimization both in the pupil plane and the image plane. In 
Section |4] we display various optimized circular pupil masks with a central obscura- 
tion and spiders. For example, we consider the pupil geometries of the 8m Subaru 
Telescope and the 3m telescope for the planned Japanese Space Infrared Telescope for 
Cosmology and Astrophysics (SPICA). Section [5] shows that we can also optimize for 
a segmented mirror; we use the pupil of the James Webb Space Telescope (JWST) 
as an example for such masks. Section |6] discusses the possibilities offered by these 
optimizations, as well as their limitations. 

2 Formalism of the two dimensional pupil optimization 

What makes these optimizations possible is a fast and accurate computation of the 
Fourier Transform (FT) of the pupil plane apodization. Indeed, the amplitude distribu- 
tion in the image plane is proportional to the Fourier transform of the electric field in the 
pupil plane, which is the arriving electric field (to be of constant, say unit, amplitude) 
multiplied by the apodization function A{x, y): 

2i7rf rD/2 .D/2 

E{u,v) = —- / A{x,y)e-^'^^dxdy. (1) 

T'^^ J-D/2J-D/2 
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Here, the aperture has diameter D along the x and y axes (the aperture is therefore 
assumed to be contained within a D x D square), the focal length of the instrument is 
F, and A denotes the light's wavelength. 

An important first step in many optical design problems is to replace quantities 
representing physical lengths with quantities that are unitless scalars times some natural 
physical scale for the problem. For example, in the pupil plane we replace x and y with 
new variables, for which we use the same notation, that measure length in units of D 
(formally, x — xD and y — >■ yD). That is, x = 1 means one diameter Similarly, 
in the image plane we replace the physical lengths u and v with values that represent 
radians on the sky (formally, u — > uXF/D and v — > vXF/D). With this substitution, 
an image plane variable u represents a physical position of u\F/D on the detector. It 
also represents an angular position of uX/D radians on the sky. Throughout the rest 
of this paper we view image plane variables as measures of \/D radians on the sky. 
Making these substitutions, ([T} becomes 

r)2 2i-^!^ 1-1/2 ,1/2 

E{u,v)^—— / A{x,y)e'^'<-^+y-Uxdy. (2) 

'^^^ J -1/2 J -1/2 

This two-dimensional Fourier transform can be approximated using discrete sets of 
points in both the pupil and the image planes: 

r)2p2»7rf Ny 

* 1=1 j=i 

with k indexed on the set {1, ... , Mu} and I on the set {1, ... , My}. We will assume 
equal sizes along the x and y axes, as well as along the u and v axes, making A^^; = 
Ny = N and Mu = My — M. Furthermore, we will consider a uniform partition in 
the {x, y)-plane, and thus Axi — Ay^ — l/N\/i. 

For the previous one-dimensional problems (where only single Fourier transforms 
were considered with the result extended to two dimensions by symmetry), the com- 
plexity of the problem was NM, where N and M are respectively the number of 
"pixels" in the pupil and image plane. A basic 2D FT leads to a complexity of N'^M'^; 
for a typical desktop computer, the maximum value for NM is roughly 4500, which 
puts serious limits on how large A^ and M can be. For example, if M is chosen to be 
30, which is about the minimum one can get away with in the image plane, then A^ can 
be at most 150. Consequently, the smallest detail in a pupil mask of 1 cm in diameter 
would have a width of ss 33/im. This is too low a resolution to produce sufficient 
contrast, as explained in more detail in Section[3] 

There is, however, a way to increase the resolution with which we sample the pupil 
plane. As described in | \1\ and 1 18], the complexity can be reduced to N'^M + NM"^ 
if the 2D Fourier transform is made in two steps, first along one axis, 

N 

Ete^p{uk,yj) = ^v4(a;„y,)e-2^'^^'"''Ax, (4) 

i=l 
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and then along the other one, using as inputs the resuhs of the first computation in the 
second ID FT, 



Eiuk,vi) = .^^ ^e-^"^^""£;tempK,^j-)Ay. (5) 
i=i 

This two-step method dramatically improves the complexity problem by reducing the 
number of computations at the expense of increased memory allocation (the interme- 
diate matrix Etemp has to be stored during this process). 

Since most telescope pupils have one or two axes of symmetry, we can use this 
feature to work with arrays 2 or 4 times smaller Additionally, because the binary 
apodization is real, the exponential term in the Fourier transform can be replaced with 
a cosine. When computing the FT this way, N can be as high as 1000, while M can be 
set to 50, and the smallest detail in the same 1 cm pupil would now have a length. 
Note that, in this particular case, we gain a factor 1000^ x 50^/(1000^ x 50 + 1000 x 
50^) « 48 in efficiency over the classical, brute-force, computation of the 2D FT. 

The maximization of the pupil mask's throughput, subject to contrast constraints, 
can be formulated as a linear programming problem: maximize E(0,0) subject to: 

-10-/2 < < 10-/2 

- E{0,0) - (6) 

< A{x^,yJ) < 1 

where 10^*^ is the targeted contrast. The electric field in the image, E{uk,vi), is subject 
to this constraint only in a specific discovery space (defined in particular by an IWA 
and an OWA). The apodization, A{xi,yj), is also defined on a specific domain, the 
initial aperture shape. 

In practice, the optimization problem is expressed in the AMPL language 1 15] and 
solved with R. Vanderbei's LOQO solver lfT6l . Figure [T] shows the AMPL model for 
the mask shown in Figure |2] 



3 Circular unobstructed aperture 

To begin, we consider a circular unobstructed aperture. Almost all coronagraphs being 
considered for a space mission require unobstructed apertures to achieve high contrast, 
and therefore necessitate complicated off-axis telescope designs. Starting with such a 
geometry allows us to compare the new approach to previously-designed shaped pupils. 
In later sections, we explore designs for more interesting pupil geometries. 



3.1 Case of a circular symmetric high contrast region 

In Vanderbei et al.f9l we found optimal shaped pupils for open circular apertures con- 
sisting of concentric rings by solving for the ring width and spacing along a radius. This 
allowed the problem to be solved in only one dimension. It is, however, interesting to 
compare the result of the ID optimization to the new 2D approach. In the latter case. 
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param pi := 4*atan(l); 

param rhoO := 3; # inner working angle 
param rhol := 15; # outer working angle 
param c := 6; # contrast 

param N := 500; # pupil plane discretization parameter 
param dx := 1/ (2*N) ; 
param dy := dx; 

set Xs := setof {j in 0.5..N-0.5 by 1} j/ (2*N) ; 
set Ys := setof {j in 0.5..N-0.5 by 1} j/ (2*N) ; 
set Pupil := setof {x in Xs, y in Ys : x"2+y"2 < 0.25} (x,y); 

var A {X in Xs, y in Ys : x"2 + y"2 < 0.25} >= 0, <= 1; 

param M := 50; # image plane discretization parameter 
set Us := setof {j in 0..M} j*rhol/M; 
set Vs := setof {j in 0..M} j*rhol/M; 
set DarkHole := setof {u in Us, v in Vs : 

u"2+v"2>=rho0'2 
&& u" 2+v" 2<=rhol ' 2 } (u,v); 

var C {u in Us, y in Ys}; # this is E_temp 

var E {u in Us, v in Vs } ; 

var area = sum { (x,y) in Pupil} A [x, y ] *dx*dy; 

maximize throughput: area; 

subject to C_def {u in Us, y in Ys}: 

C[u,y] = 2*sum {x in Xs : (x,y) in Pupil} 

A[x,y] *cos (2*pi*x*u) *dx; 

subject to E_def {u in Us, v in Vs } : 
E[u,v] = 2* sum {y in Ys} 

C[u,y] *cos (2*pi*y*v) *dy; 

subject to sidelobe_pos {(u,v) in DarkHole}: E[u,v] <= 10 " (-c/2 ) *E [ , ] ; 
subject to sidelobe_neg {(u,v) in DarkHole}: -10 " (-c/2 ) *E [ , ] <= E[u,v]; 

solve ; 

printf {(x,y) in Pupil}: "%10f %10f %10f \n", x, y, A[x,y] > "A.out"; 
Figure 1: AMPL model for the mask shown in Figure [2] with = 1000. 
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Figure 2: Pupil mask with a transmission of 36%. Top: N=200. Bottom: N=1000 (in 
both cases M=100). Left: optimal pupil mask for an IWA of 3 A/£> and an OWA of 15 
X/D. Right: magnified region contained in the red dashed rectangle displayed on the 
left. 

the optimization produces a "wiggly" version of a concentric ring mask (the wiggles 
result from the rather coarse discretization of the image plane). 

Figure |2] displays the results of the optimization with N=200 and with N=1000 
(and M=100 in both cases). In fact the optimization uses arrays 4 times smaller, taking 
advantage of symmetry, but the numbers of points we indicate afterwards are always 
the total number of points. In both cases the contrast is 10^^ inside a ring with an IWA 
of 3 \/ D and an OWA of 15 A/I?, as can be seen in Figure |3] A comparison of the 
same region is shown; the size of the arrays has a clear influence on the quality of the 
optimal transmission. In particular, the value of the transmission tends to converge to 
either or 1. The number of pixels for which the transmission adopts a value between 
and 1 decreases greatly when the size of the array is increased; the low resolution 
pupil has 1.8% of its transmitted intensity between 0.1 and 0.9, while less than 0.2% 
of the high resolution pupil has this property. 

Figure [3] shows a comparison between the PSFs of the pupils displayed in Figure |2] 
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Binary masks PSF - N=200 Binary masks PSF - N=l 000 
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Angular distance - A/D Angular distance - A/D 

Figure 3: Cross section of the PSFs of the masks displayed on Figure [2] (left: N=200, 
right: N=1000) before and after the transmission values have been rounded. The PSF 
of the optimal mask is represented by a dashed line while a red solid line is used for 
the binary version of the same mask. In both cases the intensity is represented on a 
logarithmic scale. 



with and without the transmission values rounded. In the former case, the pupil be- 
comes a true binary mask and some unwanted features brighter than 10~^ can be seen 
inside the high contrast ring when the pupil is too poorly sampled (N=200). The same 
does not occur with the higher resolution (N=1000). All of the PSFs that we show in 
the rest of the paper are computed for the binary masks obtained by rounding all pixels 
either to zero or to one. 

There is, however, one exception. We do the comparison between a true concen- 
tric ring mask (illustrated in Figure 4 of f9\) and a similar mask (see Figure |4]i fully 
optimized in 2D for the same IWA and OWA (3.5 and 7 A/D) and the same contrast 
(IQ-iOy -pjjg original concentric ring mask and the new one have similar throughputs 
(23.5% instead of 25%) but the number of rings is not the same in the two cases (5 
instead of 6). As seen in Figure |5] the mask (designed with N=1000, M=100), only 
provides a mean contrast of 2 x 10^^ when its transmission values have been rounded. 
Less than 0.12% of the pixels have an apodization value between 0.1 and 0.9, but the 
targeted high contrast level makes this small ratio significant enough to change the 
effective contrast by more than an order of magnitude. We should, however, point 
out that the ID mask was designed using a two-step process. In the first step, a lin- 
ear programming problem was solved to find a discretized approximate solution to the 
problem. For the second step, the on-off and off-on thresholds from the approximate 
solution were used as input to a ID nonconvex nonlinear optimization problem that 
hones these thresholds to highly precise values. We have not adapted this second stage 
to the current 2D design process. We are currently working on such an enhancement. 

It is also noteworthy that the new pupil masks fail to present the same circular 
symmetry as in the ID case, although they share strong similarities. This is explained 
by the fact that the PSF is computed with a small number of points, M=100. Figure |6] 
illustrates how this effect is reinforced if M is set to 50. 

The transmission of the mask does not significantly change with different values of 
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Figure 4: Left: Pupil mask designed for a contrast of 10^^". Optimization parameters 
are N=1000, M=100. The throughput is 25%. Center and right: Detail of the mask 
before and after rounding the transmission. 




8 




Figure 6: Two pupil masks optimized with a different sampling precision of the image 
plane (N=1000 in both cases while M=100 on the left, and 50 on the right). Constraints 
on the contrast and the IWA and OWA are the same as for the pupils in Figure |2] 

M and N, as long as the IWA and OWA remain the same. Therefore, it is possible to 
compute a first mask with a low resolution (in about 15 minutes for N=200, M=50). If 
the first estimate of the throughput is deemed high enough, a second optimization can 
be launched for a much larger array (which will take several hours for N=1000, M=50 
on a modern desktop computer). 

3.2 Case of two symmetric dark holes 

Instead of an annular dark hole, one may wish to optimize the pupil mask so as to 
provide high contrast in more spatially-Umited regions of the image plane. This alter- 
native is particularly relevant when the mask is to be used with a dedicated wavefront 
correction system that can only create small dark holes with high contrast. In particu- 
lar, as described in ||19l . masks can be designed with dark holes at a lower contrast and 
then paired with two DMs to produce the remaining contrast. This pairing results in an 
overall system with higher throughput and smaller inner working angle. 

Several tradeoffs exist between the IWA, the throughput and the contrast. Two 
are illustrated in Figure |7] A comparison between the pupil displayed on the first and 
second rows indicates that a smaller IWA can be obtained when the size of the dark hole 
is reduced. The reader will note that, while the throughput remains almost the same, 
the overall shape of the mask changes dramatically. Furthermore, if a higher contrast 
is sought, the throughput may be conserved if the IWA is increased, as is illustrated by 
the pupils on the first and the third rows. Although not illustrated here, reducing the 
area of the dark hole may also help increase throughput. 
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Figure 7: Pupil masks (left) design to reach 10~^ (first two rows) and lO"'^ (third row) 
with a circular unobstructed aperture, and thek associated PSFs (right). In both cases 
the high contrast region is a ring section with an OWA of 15 \/D. The IWA are 2.5 
X/D (top), 3 X/D (middle) and 4 X/D (bottom). The dark holes have an angular 
extension of 30 degrees in the first case and 45 degrees in the last two. The PSF are 
computed for true binary masks, the original transmissions being artificially rounded. 
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Figure 8: Pupil mask (left) designed to reach 10 ^ with the SPICA telescope. The high 
contrast areas can be seen on the associated PSF (right). The IWA is 3.5 \/D and the 
OWA is 12 X/D. The vertical width of the dark hole is 50 X/D. 

4 Circular pupil with a central obstruction and spiders 

Since this new two-dimensional approach can create masks for arbitrary pupil geome- 
tries, we can, in particular, find masks designed for telescopes with any central obstruc- 
tion or spider structure. This allows the coronagraph to maximize the amount of light 
collected by the telescope. In this section, we present a few masks for two real-world 
telescope pupils: SPICA and Subaru. These two pupils have two axes of symmetry, 
allowing us to simplify the computation as explained in section[2] 

4.1 SPICA 

SPICA, the Space Infrared Telescope for Cosmology and Astrophysics, is a planned 
infrared space telescope consisting of a 3m primary mirror and a secondary mirror of 
0.6 m. The central obstruction ratio is then 0.2. The spider's thickness is set to 0.08 m. 
It should be noted that the design of SPICA is subject to change and the primary and 
secondary mirrors could respectively have a diameter of 3.2 m and 0.8 m. In the latter 
case, this would have important consequences on the performances of the coronagraph 
in terms of throughput, IWA and contrast. 

Several masks for SPICA, as well as for Subaru and JWST, have been described in 
II20I I2TI. Different strategies have been explored that favor the contrast over the IWA or 
vice-versa. The mask we present here, shown in Figure[8] has a throughput of 55%. Its 
IWA is 3.3 X/D and its OWA is 12 X/D. The discovery zone for which it is designed is 
similar to the one created by the first mask that appears in \2\\. Because of the absence 
of any constraints on the optimization of the mask, a significant gain in throughput is 
obtained (the light is transmitted in an additional one fourth of the pupil surface). 
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Figure 9: Pupil masks (left) designed to reach 10^^ (top) and 10~^ (bottom) with the 
Subaru telescope. The high contrast areas can be seen on their associated PSFs (right). 
The IWA are 3 and 5 A/Z?, the OWA are 16 and 20 A/Z3, and the vertical widths of the 
dark holes are 64 and 60 X/D. 

4.2 Subaru 

UnUke SPICA, the spiders of Subaru are not perpendicular. A solution proposed in 
II20I was to place two barcode masks in a direction perpendicular to one of the spiders 
so that the design of the mask can be affected as little as possible. This way the mask 
can provide a contrast of 10^^ at 3 X/D for a throughput of 24%. 

The first mask we present here was designed with the same constraints as Mask-1 
in II2OI but rotated by w 38 degree. The particular orientation chosen in their case is 
necessary so as to display the barcode masks horizontally. We are not subject to the 
same requirement, and we have chosen to align the axes of symmetry of the central ob- 
struction with the X and y-axes of the simulation window. The top of Figure|9]displays 
the mask and its PSF It allows through almost twice the light, with a throughput of 
45%. 

The second mask was designed to reach a higher contrast, 10^*, although a larger 
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Figure 10: Pupil mask (left) designed to reach 10 ^ with the JWST. The hexagonal 
high contrast area can be seen on its associated PSF (right). The IWA and OWA are 5 
and 15 \/D. 



IWA of 5 X/D is necessary for the throughput to remain comparable (35% of the light 
is transmitted). 



5 Asymmetric pupil with a segmented mirror 

We now apply the mask optimization to the James Webb Space Telescope (JWST) 
pupil. The pupil mask is optimized with N=1024. Contrary to all the preceding pupils 
we have considered, this pupil only has one axis of symmetry. Because of that we 
cannot optimize it in the exact same way as we did for the others. We still take ad- 
vantage of one axis of symmetry in the computation of the first Fourier transform (and 
one should make sure that the pupil has the correct orientation, or rotate it by 7r/2 if 
necessary), but the second FT must be completed without it. This results in an array 
twice as large as with the other pupils (512 x 1024). 

Another consequence of the asymmetry is that the exponential term in the expres- 
sion of the FT cannot be changed into a cosine. Hence each number has a real and 
an imaginary part, which doubles the amount of memory required to solve this model. 
Furthermore, to keep the problem linear, the natural contrast constraint \E{u,v)\ < 
10~^-^E{0, 0) is replaced with a pair of constraints, one for the real part and one for 
the imaginary part: 

-lO-^-^E{0,0)/V2 < ^{E){u,v) < lO-^-'^E{0,0)/V2 
-10-^-^E{0,0)/V2 < ^{E){u,v) < 10"^-^i;(0,0)/^. 

The JWST pupil is particularly interesting: the internal and external hexagonal 
segments making up the pupil have a great effect on the contrast that can be reached. 
A first mask was generated for a circular ring dark hole with a minimal IWA of 5 X/D. 
A smaller IWA (4 X/D) could be obtained with a similar throughput by changing the 
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dark hole to a hexagonal ring, with internal and external contours rotated by 7r/2 with 
respect to each other. The resulting mask and PSF is shown in Figure 10 Several other 
geometric aspects still remain to be tested. 



6 Final Remarks 

Previous two-dimensional shaped-pupil optimizations all assumed specific geometries 
in order to reformulate the 2D problem as a ID mathematical description that could 
then be solved. We have shown in this paper that full optimizations of pupil masks in 
two dimensions are possible. No constraint other than the initial geometry of the aper- 
ture are placed on the pupil. This allows us to find new, and more optimal, solutions 
to a variety of high contrast problems. We have in particular compared composite bar- 
code masks designed for the SPICA coronagraph and for the Subaru telescope to new, 
fully optimized 2D masks. Given the same discovery zones (same contrasts, IWA and 
OWA), higher throughputs can always be expected. More importantly, these through- 
puts are simply the highest possible. 

Complex pupil geometries can be treated without much more efforts than basic 
geometries. Masks optimized for the JWST illustrate particularly well the possibilities 
that full 2D optimizations offer. Despite the central obstruction, the spiders and the 
inter-segments gaps, optimal mask solutions are found. 

Another very important observation is that the result of the optimization always 
tends to be a binary mask, and that, given a sufficiently high enough number of points 
used to sample the pupil plane, it virtually is. This had already been observed when 
optimizing ID masks, but the optimization of ripple masks and star shaped masks 
assumed a binary amplitude distribution. 

One can argue that these optimizations are constrained by the choice of the high 
contrast regions. In the basic case of a circular symmetric pupil, without central ob- 
struction and spiders, many types of dark holes can be considered (circular symmetric 
rings, portions of rings, rectangular areas, etc.), but the particularities of pupils with 
more complex features must be taken into account when designing the dark holes. For 
example, the small size of spiders leads to diffraction effects on a large scale, and the 
geometry of the spiders must be used to design the dark holes. A system of 4 spiders set 
in an orthogonal fashion will create diffraction spikes in a similar orthogonal fashion. 
One can then choose to move the diffracted energy along only one axis, so as to obtain 
a good contrast on the other, or one can also choose to maximize the contrast in the 4 
quadrants of the image plane by concentrating the diffracted energy along both axes. 

Many coronagraphaMare meant to be used with telescopes without central obstruc- 
tions. Optical devices can be used to artificially remove the spiders and/or the central 
obstruction. Masks could be designed, not to aim for high contrasts, but to attenuate 
the diffraction spikes. The advantage of choosing this solution would be a better achro- 
maticity and little to no additional aberrations introduced in the wavefront. This will 
be explored in a future paper 

The masks that are presented here are all binary, but none of them is fully struc- 
turally connected, in contrast to bar code masks, star shaped masks, or shaped pupils. 

' in particular the PIAA, the FQPM, and the vortex coronagraph 
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This makes it mandatory for them to be laid on glass or mirrors. Additional ampli- 
tude and phase aberrations will thus be introduced on the wavefront, and will have to 
be taken into account in the design of the wavefront correction system. Checkerboard 
masks are also not structurally connected, and some were manufactured for the SPICA 
project using e-beam lithography of aluminum coated BK7 glass substrates. Masks de- 
signed for 10~^° were tested in the laboratory, and contrasts smaller than 7 x 10~^ were 
achieved in laser light [21 1. More recent experiments were conducted in broadband, 
and in a vacuum chamber (see 1221 ). This time the contrast was limited by ghosts 
and ranged from 3.5 x lO"'^ for a central wavelength Aq = 650 nm (AA = 8 nm) to 
2.6 X 10~^ for Aq = 850 nm (AA = 55 nm). In both experiments the wavefront was not 
actively corrected. 

We intend in the near future to have some masks manufactured in a similar fashion, 
and use them in our high contrast imaging laboratory in Princeton, together with our 2 
DM wavefront correction system. 
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